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The  refractive  index  structure  function  coefficient,  C  ,  is  an 

n 

atmospheric  parameter  needed  to  describe  scintillation  and  small-scale 
phase  fluctuations  of  electromagnetic  radiation  propagated  in  the 
atmospheric  surface  layer.  Since  systematic  direct  measurements  of  Cfi 
for  many  climates  and  seasons  are  not  available,  an  indirect  method  is 
developed  where  C^  is  calculated  from  the  estimates  of  sensible  and 
latent  heat  flux  components  of  the  surface  energy  budgets.  This  indirect 
method  is  primarily  for  heights  less  than  4  meters,  because  low  inter - 
mittency  and  I^/K^  =  1  are  assumed.  Diurnal  variations  of  Cn  at  several 
heights  above  land  for  six  combinations  of  climates,  seasons,  and  surface 
conditions  are  calculated  from  heat  fluxes  measured  by  different  investi¬ 
gators  at  many  locations,  for  moderate  to  high  wind  velocities.  These 
predictions  of  Cn  agree  well  with  some  direct  measurements  of  Cn  when 
assumptions  of  nearly  ideal  weather  and  sites  are  met.  The  effect  of 
water  vapor  on  is  usually  a  reduction  of  about  3/3  per  cent  and  thus 
is  usually  negligible  above  land,  but  can  be  significant  above  tropical 


oceans . 
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INTRODUCTION 


For  optical  radiation  the  atmospheric  index  of  refraction  variations 

primarily  result  from  temperature-induced  density  changes.  These  variations 

cause  scintillation  and  phase  fluctuations'  of  laser  beams  propagated  in 

the  atmosphere  and. the  loss  of  resolution  of  objects  viewed  through  high- 

magnification  telescopes.  Temperature  variations  also  help  cause  the 

scattering  of  acoustic  and  radar  signals,  phenomena  used  in  the  remote 

sensing  of  the  atmosphere.  The  length  scales  of  the  temperature 

inhomogeneities  that  distort  image  and  laser  beam  propagation  in  the 

atmospheric  surface  layer  are  usually  on  the  order  of  centimeters  for 

scintillation  and  tens  of  centimeters  for  phase  fluctuations,  and  the 

eddies  with  these  sizes  are  in  the  inertial  subrange  where  fluctuations 

are  isotropic.  Isotropy  greatly  aids  the  development  of  theories  like 

those  in  Tataraki  [1961]  that  describe  electromagnetic  propagation,  and 

only  one  atmospheric  parameter,  the  refractive  index  structure  function 
2 

coefficient,  Cn,  is  needed.  Using  Kolmogorov's  theory  one  can 

determine  for  optical  radiation  from  a  variance  of  temperature  \ 

differences  between  two  fast-response  thermometers  or  from  a  temperature 

power  spectrum  at  one  thermometer.  When  the  ideal  case  of  the 

temperature  and  wind  structure  being  horizontally  homogeneous  in  the 

surface  layer  exists,  one  measurement  suffices  and  Cn  can  be  inferred 

for  any  height  within  a  few  meters  of  the  surface. 

The  purpose  of  this  paper  is  to  present  diurnal  variations  of 
for  optical  radiation  at  several  heights  in  the  atmospheric  surface 
layer  for  various  climates,  seasons,  and  surface  characteristics. 

However,  since  direct  systematic  measurements  of  Cn  are  not  available, 
an  indirect  method  is  developed,  compared  with  similar  indirect 
methods,  and  used  to  obtain  diurnal  trends.  The  indirect  method  uses 
sensible  and  latent  heat  flux  components  of  the  surface  energy  budgets, 
which  have  become  abundantly  available  in  the  meteorological  and 
climatological  literature.  Representative  wind  velocities  are  chosen 
and  the  effect  of  water  vapor  is  evaluated.  Also,  the  model  developed 
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here  is  for  use  primarily  at  heights  less  than  4  meters,  since  the 
effect  of  intermittency  at  greater  height  limits  the  validity  of 
necessary  empirical  functions. 

Some  of  the  terminology  used  in  this  paper  must  be  explained. 

First,  "turbulence"  is  not  synonymous  with  "temperature  fluctuations" 
as  is  common  in  atmospheric  optics  [Lawrence  et  al*»  1970].  It  will  be 
shown  that  water  vapor  fluctuations  can  be  important  relative  to 
temperature  fluctuations  above  a  lake  or  an  ocean.  Also,  wind  velocity 
fluctuations  can  significantly  affect  the  propagation  of  high-power 
laser  beams.  Thus  "turbulence"  implies  any  fluctuating  property  of  the 
air,  and  equating  a  value  of  with  intensity  of  turbulence  is  mis¬ 
leading  unless  thermal  turbulence  is  specified.  Second,  "stable"  con¬ 
ditions  refer  to  an  inversion  or  an  increase  of  potential  temperature 
with  height,  and  heat  being  convected  from  the  air  to  the  cooled  ground. 
"Neutral"  conditions  are  adiabatic  conditions;  approaches  zero  and 
the  average  heat  flux  is  zero.  "Unstable"  conditions  are  lapse 
conditions,  or  a  decrease  beyond  the  normal  adiabatic  decrease  of 
temperature  with  height,  and  heat  is  being  convected  from  the  surface 
to  the  air. 


INDIRECT  CALCULATION  OF  C 

n 


by 


The  structure  function,  D^,,  for  temperature  fluctuations  is  given 


Dt  =  <(TrT2)2>, 


(1) 


where  T^  and  T2  are  the  temperatures  at  two  points  fixed  on  a  line 

oriented  normally  to  the  mean  wind  direction  and  separated  by  distance 

r.  1\tfo  thermometers  can  be  used  to  measure  T^  and  T2,  but  an 

alternative  method  is  to  use  one  thermometer  and  determine  T.,  by 

«« 

lagging  in  time  from  T^  by  At.  Then  Taylor's  Hypothesis  of  a  frozen 
field  is  assumed  and  r  =  At  u,  where  u  is  the  mean  longitudinal  wind 
speed  [Lmley  and  Panofsky,  1964].  When  r  is  of  the  order  of  inertial 
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subrange  scales,  Kolmogorov's  reasoning  can  be  applied  [ Obukhov 3  1949; 
Corrains  1951]  to  yield  the  following: 


where  the  coefficient  is 


Dt  =  C2r2/3, 


C2  =  ANe"1/3. 


The  value  of  A  is  constant,  N  is  the  rate  of  dissipation  of  temperature 
fluctuations,  and  e  is  the  rate  of  dissipation  of  mechanical  turbulent 
energy. 

Tataraki  [1961]  suggests  that  in  a  homogeneous,  statistically 
stationary,  atmospheric  surface  layer  the  rate  of  dissipation  of  tempera¬ 
ture  fluctuations  is  equal  to  the  product  of  the  vertical  heat  flux,  H, 
and  the  mean  vertical  temperature  gradient,  3T/3z.  Data  given  by 
Wyngaard  et  al.  [1971]  indicate  that  this  is  a  good  assumption. 

N  =  «  ||  CpCp)  (4) 

where  (pC  )  is  the  heat  capacity  of  air  and  H(pC  )  =<w*T'>  is  the  co- 
P  P 

variance  of  the  vertical  wind  and  temperature  (the  primes  indicate 

deivations  from  the  means).  In  the  atmospheric  surface  layer,  e  is  the 

sum  of  mechanical  and  buoyant  production  rates  of  mechanical  and 

turbulent  energy  [Lumley  and  Panofsky ,  1964] . 

e  =  u*2  |»  U  -  Rf)  (5) 

where  u*  =  <u'w'>  '  is  the  surface  friction  velocity,  u  is  the  mean 
longitudinal  wind  speed,  and  Rf  is  the  flux  Richardson  number.  In  this 
paper,  terms  that  may  cause  equation  5  to  be  in  error  by  up  to  40  percent 
for  non-neutral  conditions  [Record  and  Cramers  1966;  Znbkovskiy  and 
Koprov,  1970;  Wyngaard  and  Cotes  1970;  MoBean  et  ai,t  1971]  are  neglected 
because  of  disagreement  about  the  sign  and  magnitude  of  the  error  terms. 
The  value  of  A  in  equation  3  is  obtained  from  temperature  spectrum 
analyses  and  a  representative  value  of  2.8  is  used  here  [ Panofsky s  1969 
Pond  et  al.s  1971;  Wyngaard  and  Cotes  1971],  This  value  is  very  close 


/ 


to  the  constant  for  the  transverse  wind  structure  function  coefficient, 
but  no  theory  states  these  constants  must  be  the  same.  Since  temperature 
is  a  scalar,  there  is  no  longitudinal  form  of  equation  2  as  exists  for 
wind  components,  where  the  longitudinal  constant  is  3/4  the  transverse 
constant. 

Substituting  equations  4  and  5  into  equation  3,  we  find 

C2  =  2.8  <w'T*>  ||  [u*2  ||  (1-Rf)  ]  1/3  .  (6) 


Lumley  and  Panofsky  [1964]  list  several  well-known  relationships  we 
can  use  to  simplify  equation  6. 


<u,"’>  =  kmI7'  (7a) 

<“'T’>  ■  "h  §  •  <7b> 

where  and  are  the  eddy  diffusivities  for  momentum  and  heat, 
respectively.  The  average  wind  speed  and  temperature  gradients  are 


3z  kz 


3T 

3z 


kz  ’ 


(8a) 

(8b) 


where  k  =  0.4  is  von  Karmon's  constant,  T*  =  (w'T’^  /u* ,  and  <j>M  and 
(fijj  are  the  nondimensional  .shears,  of  wind  and  temperature,  respectively. 
Substitutions  of  equations  7  and  8  into  equation  6  result  in  these 
final  equations : 


C^,z 


1/3 


/  T*  =  1.67$ 


1/2  iy/2,,1/3 


M 


k  (1-Rf) 


1/6  n  -1 


O) 


and 


C/tf^273)  -  1.67k2/V/2  [*J/3<l-Rf)1/6]_1,  (10) 
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where 


a  =  S/Sl  =  V^H  * 


i 

i 


j 


i 

!fe 

i 

i;. 

i 

Sft 

-  ft 


ft' 


Atmospheric  variables  can  be  nondimensionalized  with  u*  and  T*  since 
they  do  not  vary  with  height  in  the  atmospheric  surface  layer.  All 
properly  nondimensionalized  parameters  can  be  characterized  in  terms  of 
z/L  where  L  is  the  Monin-Obukhov  length. 

z/L  =  zkgT* (Tu*2)-1 ,  (11) 

where  g  is  the  acceleration  of  gravity.  The  right-hand  sides  of 
equations  9  and  10  have  three  parameters  that  must  be  evaluated  as  a 
function  of  z/L.  The  flux  Richardson  number  is  an  explicit  function  of 
z/L:  Rf  =  <{>M  z/L.  The  function  <J>M  has  been  established  by  the  empirical 
KEYPS  function,  but  for  computation  purposes  we  choose  to  use  the 
following: 


<(»M  =  1+4.5  (z/L) 

t — \ 

o 

A 

l-J 

M 

(12  a) 

4>m  =  l+4.5(z/L)exp(4.5z/L) 

(0.>  z/L  >  -0.04), 

(12b) 

<f>M  =  0.08+0.285(-z/L)1/3 

(z/L  <  -0.04). 

(12c) 

These  equations  result  in  curves  of  <f>^  verses  z/L  that  closely  agree 
with  those  given  by  Businger  et  al.  [1971].  For  a,  however,  widely 
different  estimates  have  been  obtained  [Laykhtman  and  Ponomarea,  1969; 
Businger  et  al.  t  1971].  Here  we  will  base  our  estimates  on  measure¬ 
ments  described  by  Wesely  et  al.  [1970]  which  were  within  4  meters  of 
the  surface.  The  value  of  a  was  apparently  near  unity  for  all  stability 
conditions  since  flux  estimates  from  aerodynamic  methods  using  a  =  1 
agreed  with  direct  flux  estimates  via  eddy  correlation  techniques.  The 
following  forms  for  a  are  assumed: 

a  =  (1+2.25Z/L)  (z/L  >  0),  (13a) 

a  =  l-0.225z/L  (z/L  £  0) .  (13b) 

13 


At  heights  greater  than  4  meters,  equation  13  is  inaccurate  since,  as 
shown  by  Haugen  and  Kaimal  [1971],  intermittency  increases  with  height 
and  calculations  by  Steams  [1971]  show  that  with  highly  intermittent 
conditions  a  becomes  much  larger  than  unity. 

After  substitution  of  equations  12  and  13  into  equation  9,  the  curves 
in  Fig.  1  can  be  calculated.  Also  shown  is  a  curve  determined  from  data 
collected  at  z  >  5.66  meters,  as  summarized  by  Wyngaard  et  at.  [1971]. 

The  two  curves  differ  because  of  the  higher  level  of  intermittency 
that  existed  at  the  greater  height,  thus  increasing  a  or,  alternatively, 
increasing  T*  with  respect  to  C,j,  in  equation  9. 

The  value  of  C,p  can  also  be  determined  from  the  mean  gradients  of 

temperature  and  horizontal  wind  rather  than  the  fluxes  of  heat  and 

momentum.  This  method  uses  the  gradient  Richardson  number,  Rg,  as  the 

stability  parameter:  Rg  =  ^”^a"^z/L.  Substituting  equations  12  and 

13  into  equation  10  results  in  Fig.  2,  which  also  shows  results  from 

the  same  experiment  cited  for  Fig.  1  (z/L  >_  5.66  meters).  The  effect 

of  the  increased  level  of  intermittency  at  z  >  4  meters  is  to  increase 

3T 

the  ordinate  values  because  a  is  increased  or,  alternatively,  is 
decreased  with  respect  to  Cp  in  equation  10.  The  remaining  curve  in 
Fig.  2  was  determined  by  Tsvang  [1960]  from  spectral  analysis  of  tempera¬ 
ture  fluctuation  measured  at  heights  of  1  and  4  meters.  Ts vang's  curve 
appears  low  for  unstable  conditions  (Rg  <  0) ,  as  is  also  indicated 
by  his  comparisons  with  previous  Russian  estimates  obtained  with  equation 
1.  Experimental  data  by  Livingston  et  al.  [1970]  have  shown  that 
Tsvang1 s  curve  is  reasonably  accurate. 

If  only  temperature  changes  significantly  affect  the  air  density 
and  thus  its  refractive  index, 


where  p  is  the  static  atmospheric  pressure  and  A^  is  a  constant  for  a 
particular  radiation  wavelength.  However,  fluctuations  in  water  vapor 
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Figure  1  -  Dependence  of  Cj  on  z/L.  The  solid  line  was 
computed  from  equation  9. 


Figure  2  -  Cj  nondimensionalized  with  the  local  mean 

temperature  gradient  plotted  against  Rg.  The 
solid  line  was  computed  from  equation  10. 


concentration  can  also  affect  air  density.  To  evaluate  this  effect,  we 
assume  that  the  eddy  diffusivities  for  scalar  quantities  are  identical 
at  any  given  time  and  location.  This  is  justified  by  the  success 
[e.g.,  Wesely  et  al.s  1970} > of  equating  ^  with  K^,  the  eddy  diffusivity 
for  water  vapor,  in  the  Bowen's  ratio  method  of  computing  vertical 
sensible  and  latent  heat  flux  in  the  atmospheric  surface  layer  [Lumley 
and  Panof8kyt  1964],  With  replacement  of  T  by  n  in  equations  7b  and  6 
and  the  assumption  of  Kn  =  K^,  it  can  be  shown  that 

Cn  =  CT  <w’n’>  /  <W'T'>  (15) 

The  atmospheric  refractive  index  is  given  by  Fleagle  [1950]  as 
n  -  1  =  A^p  /  T  +  k^e/1 , 

where  e  is  the  water  vapor  pressure  and  is  a  constant  for  a  particular 

radiation  wavelength.  Both  k^  and  A2  can  be  evaluated. with  the  Barrell 

and  Sears  formula  [List,  1958]  using  the  technique  described  by 

Kallistratova  and  Timanovaky  [1971];  A  =  78.7  x  10_6oK  mb-1  and 
-6  -1  ^ 

k^  -  12.4  x  10  °K  mb"  for  a  wavelength  of  0.6y.  The  value  of  A ^ 
varies  only  a  few  percent  over  the  visible  and  near  infrared  spectrum, 
but  A2  varies  greatly  for  different  infrared  wavelengths  because  of 
anomalous  dispersion  [ Chalfee ,  1968]. 

When  pressure  fluctuations  and  third  order  correlations  are 
neglected,  it  can  be  showh  that  for  visible  light 

<w'n’>  «  -A^T  "2  <w'T'>  [l-TA2Cp(0.622A1VB)"1]  .  (16) 

Therefore,  at  T  =  288°K  and  p  -  1013  mb, 

Cn  =  C^pT  "2  [1  -  0.02988"1] ,  (17) 
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where  V  is  the  latent  heat  of  vaporization  and  6  is  Bowen's  ratio. 

Since  B  =  H/E  where  E  is  the  latent  heat  flux,  a  high  value  of  B  in  the 
daytime  usually  means  that  energy  used  in  evaporation  is  much  smaller 
than  the  thermal  convective  loss  from  the  heated  surface.  This  is  the 
most  frequent  case  above  land  since  B  >  0.2  is  typical,  and  the  correction 
in  equation  17  is  negligible  there.  Usually  B  >  0  because  the  earth's 
surfaces  are  most  often  either  warming  and  evaporating  or  cooling  and 
condensing;  air  parcels  warmer  than  the  surrounding  air  usually  have  a 
higher  specific  humidity.  However,  when  surface  cooling  and  evaporation 
occur  simultaneously,  6  <  0  can  be  found  above  oceans  and  lakes  when 
severe  advection  occurs  and  occasionally  in  a  dry  atmosphere  above  moist 
ground  at  night. 


DIURNAL  CYCLES  OF  C 

n 

Diurnal  trends  of  Cfi  are  shown  in  Figs.  3  through  8  for  different 
micrometeorological  conditions  above  land  at  middle  latitudes.  Hundreds 
of  typical  trends  can  occur,  but  a  particular  micrometeorological  situ¬ 
ation  on  a  clear  and  cloudless  day  should  yield  either  one  of  the  diurnal 
trends  shown  here  or  something  between  two  of  the  graphs.  The  variables 
that  determine  which  curve  to  use  are  listed  in  Table  1,  and  they  can  be 
summarized  as  follows:  (1)  season  and  latitude,  which  determine  the 
diurnal  solar  radiation  cycle;  (2)  the  amount  of  water  available  for 
evapotranspiration  from  the  surface,  whose  main  effect  is  to  reduce  the 
surface  temperature  rise;  and  (3)  the  radiation  absorption  character¬ 
istics  of  the  surface,  usually  determined  by  the  type  and  condition  of 
vegetation,  which  is  strongly  dependent  upon  the  above  two  factors  and 
man's  modifications.  Another  factor  that  must  be  taken  into  account  is 
that  Cn  is  approximately  proportional  to  1/u*.  In  Figs.  3  through  8, 
moderate  wind  speeds  of  about  3  m/sec  at  a  2-meter  height  (u*  =  25  m/sec) 

are  assumed  when  H  <  125  W  m  ,  and  strong  wind  speeds  of  about  5  m/sec 

_2 

at  a  2-meter  height  (u*  =  40  m/sec)  are  assumed  when  H  >  125  W  m  . 
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TABLE  1.  Assumed  climatic  conditions 


Possible  Climates 


Case 

Rainfall*2 

Latitude 

Example 

Season 

1 

>25 

30°-50° 

Treeless 

area 

Winter 

2 

>25 

30°-50° 

Treeless 

area 

Winter 

3 

75-125 

30c-50° 

Wet 

grasslands 

Summer 

4 

50-100 

30°-5G° 

Plains, 

steppes 

Summer 

5 

<25 

O 

O 

t 

o 

O 

H 

Desert 

Winter 

6 

<25 

10°-40° 

Desert 

Summer 

acm/year 

b. 

frozen 


C  estimates  in  Figs.  3-8. 
n  & 


Soil  Moisture 

Subsurface  Surface  Surface 


b 

b 

Wet 

Wet 

Snow 

Wet& 

Wet^3 

Short  brown 
grass 

Wet 

Wet 

Short  green 
grass 

Wet 

Dry 

Short  green 
grass 

Dry 

Dry 

Bare  sand, 
sparse  or 

Dry 

Dry 

no  plants 

to  bottom)  for 


CASE 


meters  (top  to  bottom)  for  case 


Figure  5  -  C  vs  mean  solar  time  at  0.5, 


CASE 


meters  (top  to  bottom)  for  case 


CASE 


meters  (top  to  bottom)  for  case 


EMP  RANGE 


Winter  trends  of  Cn  are  shown  in  Figs.  3  and  4,  for  which  all  air 

and  surface  temperatures  are  assumed  below  freezing.  The  sensible  heat 

fluxes  and  thus  C^'s  would  be  different  above  a  thawing  surface.  Since 

the  values  of  H  in  Figs.  3  and  4  only  partially  come  from  direct  flux 

measurements,  the  values  of  C  are  the  least  reliable  of  the  six  cases 

presented  here.  Nighttime  H  over  snow  was  estimated  using  Neiderdorfer's 

-2 

[ Geiger 1961]  nocturnal  mean  of  about  17  W  m  .  Also,  temperature  lapse 
rates  measured  by  Geiger  [1961]  above  snow  indicate  that  -H  and  thus  Cn 
gradually  decrease  during  a  clear  night.  Since  during  the  daytime  only  a 
a  slight  decrease  of  temperature  with  height  above  snow  was  measured, 
small  estimates  of  H  were  used,  and  they  were  chosen  to  be  in  phase  with 
the  expected  solar  radiation.  The  magnitudes  of  H  at  night  in  Fig.  4  are 
smaller  than  for  Fig.  3  because  of  the  lessened  emitted  radiative  flux 
and  increased  soil  heat  flux  for  the  grassy  surface.  For  unstable  con¬ 
ditions  in  Fig.  4,  H  is  greater  than  in  Fig.  3  because  of  the  increased 
absorption  of  solar  radiation. 

As  shorn  in  Fig.  9,  measurements  were  taken  to  test  the  curves  in 
Fig.  4.  The  measurements  were  made  above  a  level,  brown  grass  surface 
on  the  Electromagnetic  Propagation  Range  (EMP  Range)  at  the  U.  S.  Army 
Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  Maryland.  The 
curves  in  Fig.  4  are  shown  to  be  reliable  in  Fig.  9,  despite  the  low 
values  of  Cn  between  1700  and  2400  hours  caused  by  extremely  low  winds, 
a  common  but  sometimes  unpredictable  occurrence. 

In  Fig.  5,  Cn  trends  are  shown  for  a  wet,  grassy  surface  in  late 

spring  or  late  summer,  or  in  midsummer  with  a  haze.  The  values  of  H 

chosen  are  typical  of  Central  Europe  [Geigerj  1961]  and  were  obtained 

from  measurements  taken  near  Hancock,  Wisconsin, by  Wesely  et  al.  [1970]. 

In  Fig.  10  measurements  of  Cn  above  grass  at  the  EMP  Range,  Aberdeen 

Proving  Ground,  Maryland, show  that  the  predictions  are  accurate.  Clouds 

caused  low  C  's  from  1600  to  1800  hours.  Predictions  of  C  for  a  grassy, 
n  n 

partially  dry  surface  are  shown  in  Fig.  6,  where  the  values  for  H  were 
chosen  from  Wesely's  measurements  near  Davis,  California.  The  H  and 


/ 

/ 
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thus  estimates  are  typical  for  the  Great  Plains  of  the  United  States 
[Wyngaard  et  al.3  1971;  Lettau  and  Davidson 3  1971]  or  the  Russian 
Steppes  [Kallistratova  and  Timanovskiy 3  1971;  Zilitinkevich  and  Chalikov3 
1968]  provided  the  soil  is  wet  beneath  the  top  few  centimeters  and  the 
grass  is  actively  growing.  The  maximum  nighttime  C^'s  for  both  Figs.  5 
and  6  occur  an  hour  or  two  after  sunset  because  evaporation  continues  to 
cool  the  surface. 

Figs.  7  and  8  are  extremely  dry  climates  in  the  winter  and  summer, 

respectively.  However,  the  values  of  H  and  Cn  above  grassy  plains  or 

steppes  during  a  summer  drought  could  approach  those  in  Fig.  7  as  the 

soil  dries  out  and  the  grass  turns  brown.  The  values  of  H  in  Fig.  7  are 

taken  directly  from  measurements  made  in  Pampa  de  La  Joya,  Peru,  by 

Steams  [1969].  Daytime  values  of  H  in  Fig.  8  were  taken  from  Vehrencamp' s 

[1953]  measurements  at  El  Mirage,  California,  but  since  his  nighttime 

values  of  wind  speed  and  H  were  low,  the  nighttime  values  of  H  and  Cn  in 

Fig.  8  were  taken  from  Fig.  7.  The  predicted  C^'s  of  Fig.  8  have  been 

partially  substantiated  by  measurements  made  by  T.  H.  Pries  at  the  U.  S. 

Army  Atmospheric  Sciences  Laboratory  at  White  Sands  Missile  Range  where 

preliminary  evaluation  of  data  collected  in  early  1972  show  peaks  near 
-7  -1/3 

C  =  7  x  10  m  at  a  3-meter  height, 
n 

The  trends  drawn  in  Figs.  3  through  8  are  poor  climatological  models 
because  several  assumptions  of  "ideal"  conditions  are  made.  Large  re¬ 
ductions  in  C  occur  whenever  clouds  obscure  the  sun  or  a  significant 
n 

portion  of  the  sky  at  night.  Reductions  in  visibility  also  lower  the 

radiation  absorbed  at  the'  surface  and  thus  the  C  ..  Real  diurnal  curves 

n 

are  rarely  as  smooth  as  these  "ideal"  ensemble  averages  since  wiggles, 
usually  caused  by  nonstationary  winds  and  other  synoptic  effects,  can 
be  expected.  Equations  9  and  10  are  valid  only  when  the  site  is  flat 
and  homogeneous  upwind  from  the  measurement  site.  Estimates  of  the 
necessary  fetch  have  been  determined  by  several  authors  [Elliot3  1958; 
Brooks,  1961;  Panofeky  and  Townsend,  1964;  Taylor ,  1969]  to  be  from  50 
to  100  times  the  height  of  the  sample.  If  intermittency  is  strong,  which 
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is  most  likely  to  occur  at  heights  greater  than  4  meters,  then  equation 
13  is  invalid  and  Cn  will  be  overestimated  by  equations  9  or  10.  Finally 
discontinuities  of  Cn  do  not  exist  as  shown  on  the  graphs  for  the 
neutral  conditions  near  sunrise  and  sunset.  The  curves  should  approach 
zero,  but  Cn  is  usually  finite  because  the  "fleeting  moment"  at  which 
there  are  no  temperature  fluctuations  is  included  in  a  ten-to-thirty- 
minute  average  essential  for  point  measurements  in  the  atmospheric 
surface  layer.  The  curves  for  different  heights  should  be  slightly  out 
of  phase  because  of  Swinbank's  [1964]  observation  that  the  moment  at 
which  the  temperature  fluctuations  go  to  zero  at  the  surface  may  be  as 
much  as  15  minutes  sooner  than  at  16  meters. 

CONCLUSION 

The  diurnal  cycle  predictions  of  C^  presented  in  this  paper  indi¬ 
cate  the  magnitude  of  refractive  index  fluctuations  in  the  inertial 
subrange  of  the  atmospheric  surface  layer  over  land.  These  values  of 
Cfi  can  be  used  in  the  various  theoretical  expressions  that  give  the 
magnitude  of  amplitude  and  phase  alterations  of  electromagnetic  radiation 
through  a  thermally  turbulent  atmosphere.  Low  intermittency  is  assumed, 
and  thus  the  results  are  most  applicable  to  heights  less  than  4  meters. 
For  a  given  climate,  season  and  surface  type  the  values  of  Cn  given  here 
are  the  maximum  expected  since  cloudless  skies,  excellent  visibility  and 
nearly  optimum  wind  are  assumed.  Thus  the  diurnal  cycles  in  Figs.  3 
through  8  should  be  used  only  as  rough  approximations  for  clear  skies, 
and  during  propagation  experiments  measurements  of  Cn  should  be  made. 

The  simple,  direct,  structure  function  approach  of  equation  1  is  more 
desirable  than  the  indirect  methods  using  equations  9  and  10,  which 
rely  upon  semiempirical  relationships.  Regardless  of  the  measurement 

technique  chosen,  one  measurement  usually  suffices  to  determine  C  at 

-1/3 

any  height,  since  as  shown  by  equation  9,  Cn  is  proportional  to  z 
in  the  atmospheric  surface  layer. 
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Future  work  in  determining  the  diurnal  behavior  of  in  the 
atmospheric  surface  should  include  a  study  of  the  influence  of  water 
vapor  fluctuations  above  oceans,  lakes,  and  swamps.  Recent  flux 
measurements  by  Pond  et  al.  [1971]  above  the  sea  show  3  0.1  within 

10°  of  the  equator,  and  according  to  equation  10  an  estimate  of  based 
upon  a  direct  measurement  of  would  be  too  large  by  30%.  A  smaller 
effect  is  expected  at  greater  latitudes  since  seasonally  averaged 
values  of  3  increase  to  nearly  0.5  at  70°N  [ Sellers ,  1965],  The  value 
of  3  above  oceans  near  shores  and  above  small  bodies  of  water  can  be 
exceptionally  small  and  is  usually  negative  during  the  daytime  if  the 
nearby  land  regions  are  dry  [ Geiger ,  1961}  Sellers,  1965].  Thus,  at 
selected  times  and  places  water  vapor  fluctuations  can  greatly  influence 
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